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ABSTRACT 

We examine colors from 3.6/im to 24/im as a function of metallicity (O/H) for a sample of 34 
galaxies. The galaxies range over 2 orders of magnitude in metallicity. They display an abrupt 
shift in the 8/im to 24/xm color between metallicities 1/3 to 1/5 of the solar value. The mean 8/Ltm 
to 24/im flux density ratio below and above 12 + log{0/H) = 8.2 is 0.08 ± 0.04 and 0.70 ± 0.53, 
respectively. We use mid-infrared colors and spectroscopy to demonstrate that the shift is primarily 
due to a decrease in the 8/im flux density as opposed to an increase in the 24/mi flux density. This 
result is most simply interpreted as due to a weakening at low metallicity of the mid-infrared emission 
bands usually attributed to PAHs (polycyclic aromatic hydrocarbons) relative to the small-grain dust 
emission. However, existing empirical spectral energy distribution models cannot account for the 
observed short- wavelength (i.e., below 8/xm) colors of the low-metallicity galaxies merely by reducing 
the strength of the PAH features; some other emission source (e.g., hot dust) is required. 

Subject headings: galaxies: ISM — infrared: galaxies 
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1. INTRODUCTION 

Early galaxies must have formed at a very low metal- 
licity. Such galaxies likely had different properties from 
typical galaxies observed in the local universe, where gen- 
erations of star formation have enriched the ISM (in- 
terstellar medium). These early galaxies are difficult to 
study because they are so distant. As a result, nearby 
low-metallicity galaxies have generated considerable in- 
terest as local analogs of young galaxies at high redshi ft 
(see, for example, the review bv lKunth fc Osthnl feOOO'). 

Understanding the infrared properties of nearby 
low metallicity galaxies will help guide interpreta- 
tion of Spitzer observations of high-redshift ones (e.g ., 
iFazio et all 120041: lYan et al.ir2004t iLagache et al.i r20041. 
An important aspect of infrared galaxy spectral en- 
ergy distributions (SEDs) is the behavior between 3 and 
24/im, a range that includes emission from PAH (Poly- 
cyclic Aromatic Hydrocarbon) molecules and warm dust. 
At redshifts between 1 and 3, the PAH features move 
through the MIPS 24um band. Ad h oc adjustments 
to th is spectral region in SED models l)Lagache et al.l 
120041) have been requi red to fit the 24 um nu mber counts 
iPapovich et al.ll2004t iMarleau et al.M2004t IcTharv et all 
120041) . iLlantTet alJ l|2004h show that infrared luminous 
galaxies at z ~ 0.6 are only ~ 50% as metal-rich as 
their local analogs. Therefore, interpreting the behavior 
of galaxies detected with Spitzer to substantially higher 
redshifts (e.g., Perez-Gonzalez et al. 2005, submitted to 
ApJ) depends critically on understanding the effects of 
metallicity on the 24/im flux densities for 1 < z < 3. 

Low-metallicity galaxies are difficult to find: I Zw 18, 
the fi rst extreme example identified l)Searle fc Sargent! 
Il972j) remains the lowest-metallicity galaxy yet found, 
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despite 3 decades of searching. Low-metallicity galaxies 
have had little or no star formation up to the present 
epoch and are therefore of very low luminosity. As a 
result, most existing compilations of the MIR proper- 
ties of galaxies are restri cted to high-metalli c ity, high- 
lumin osity targets (e.g., IGenzel et al.l 119981 ILu et alJ 
12003^) . Where low-metallicity galaxy spectra have 
been obtained, though, they differ dramatically from 
their higher-metallicity cousins in their infrared prop- 
erties, with ho tter large-grain dust emission and weak 
PAH emission jSauyage et al.lll99fi iCalzetti et al.ll2000t 



Galliano et alJ 120031: iContursi et all 1200(1 iThuan et al] 
1999t IHouck et all l2004bt IRoche et all 119911: IMadderJ 

2002). 

This work greatly expands the sample of low- 
metallicity galaxies observed in the infrared, to explore 
the effects of metallicity on the infrared properties of 
galaxies in a systematic way. The larger sample has al- 
lowed us to pinpoint the metallicity at which the transi- 
tion from weak to strong PAH emission occurs. 

2. OBSERVATIONS AND DATA REDUCTION 

The new data presented here were obtained using 
IRAC and MIPS. The IRAC data are all standard 
pipeline reductions using the most recent versions avail- 
able in the archive, ranging from 9.5.0 to 10.5.0. The 
MIPS data w ere reduced using t he MIPS DAT (Data 
Analysis Tool: lGordon et al.ir2005|) . 

Photometry was performed using the "imexam" and 
"imstat" tasks in IRAF 5 . Most sources in the sample are 
compact; in general, apertures large enough to encom- 
pass all the flux were used for MIPS, while the stan- 
dard aperture with a radius of 10 pixels (12") was used 
for IRAC, as suggested in the IRAC data handbook. 
For the extended galaxies, matching apertures that in- 
clude the whole galaxy were used for MIPS and IRAC, 
and extended-source corrections have been applied to 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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the IRAC measurements as described in the IRAC data 
handbook. No color corrections were applied, which 
should not have a strong impact on the results since the 
corrections listed in the IRAC and MIPS handbooks are 
typically only a few percent. 

To increase the number of metal-rich galaxies in the 
sample, we have supplemented our measurements with 
results from the literature. Where IRAC and MIPS mea- 
surements were available, they have been taken directly 
from the sources indicated in Table ^ We have synthe- 
sized IRAC and MIPS measurements fro m an IRS (In- 
frared Spectrogr aph: iHo uck et al] |2004a'l measurement 
of NGC 7714 by Bran dl et aLl ifeOfM) by convolving the 
spectrum with the IRAC 8/im and MIPS 24/im spectral 
response curves available on the SSC (Spitzer Science 
Center) website. 

Finally, we have added several meas urements from 
MSX (Midcourse Space Experiment) bv iKraemer et ail 
112002ft and IRAS (Infrared Astronomical Satellite) by 
iRice et alJ (|1988ft . The MSX Band A (8.28/mi) and 
IRAC 8/im flux densities measured for two galaxies 
(M 101, K. Gordon, in preparation; M 83, C. Engel- 
bracht, in preparation) are very similar (within 20% in 
both cases), so we have made no correction to the MSX 
measurements used in this paper. We have made an 
empirical color correction from the IRAS 25/im band 
to the MIPS 24/im band by comparing the m easured 
flux densities of NG C 7331 llRegan et alJl200l M 8 1 
ijGordon et al .ll2004j) . NGC 55 (jEngelbracht et al.ll200l . 
and M 101 (K. Gordon, in preparation). The measured 
24/im to 25/im ratios ranged from 0.81 to 0.94, so we 
applied the average value, 0.87, to all the IRAS mea- 
surements incorporated into the flux density ratio mea- 
surements in Table HI 

The resulting sample is a heterogeneous collection 
of star-forming galaxies without strong active nuclei. 
The global photometry measurements make this sample 
suitable for comp arison to high-redshift samples (e.g., 
iFazio et al.ll200^) . where only flux measurements inte- 
grated over the whole galaxy are feasible. The metallic- 
ity quoted for each galaxy generally applies to a much 
smaller beam than for the infrared data. The outer re- 
gions of each galaxy should have lower metallicity than 
the nuclei. Thus, the beam mismatch would, if anything, 
weaken the trend of infrared properties we see clearly in 
the suite of measurements. We conclude it is not an im- 
pediment to this study. 

3. RESULTS 

To explore the contribution of PAH emission to the 
8/im band, we compare the 8/im measurements to bands 
at both shorter and longer wavelengths. The longer- 
wavelength 24/im band is already dominated by emis- 
sion by dust and so we use it directly, but the shorter- 
wavelength 5.8 and 4.5/im bands can contain significant 
contributions from PAH emission and starlight, respec- 
tively. To avoid contamination from PAH features, we 
do not use the 5. 8/im band to derive the dust contin- 
uum measurement. We instead use the 4.5/im band, from 
which we subtract the stellar component as follows: We 
assume the 3.6/im band is dominated by starlight (which 
may result in an underesti mate of the 4.5fxm dust emis- 
sion by a few percent; cf. iHelou et alJ 12004) and mul- 
tiply the 3.6/im measurement by the 4.5/3.6/iin ratio 



expected for the stellar popu lation, derived from Star- 
burst99 ijLeitherer et aljri999f) data files available on the 
web 6 . The ratio predicted by that model is typically 0.53 
to 0.61, depending on star formation history and metal- 
licity. We adopt the average value of 0.57 and assign to it 
an uncertainty of 7%. Hereafter, we refer to this ratio as 
a. Our adopted value of a is similar, but not identical to, 
the value of .47 in the temp late elliptical galaxy spec- 
trum used by ILu et alJ <|2003) to measure the hot dust 
component in a sample of galaxies. This will have some 
effect on our measurement of the 4.5/im dust component 
for specific galaxies but does not affect the conclusions 
we draw from a diagnostic plot, as we discuss below. The 
scaled 3.6/im measurement is then subtracted from the 
4.5/im measurement. For the galaxies in this sample, the 
scaled 3.6/im measurement ranges from 26% to 96% of 
the 4.5/im measurement, with an average value of 70%. 

We did not subtract the stellar contribution from the 
8/im and 24/im measurements because the correction is 
small for these galaxies and has no impact on our conclu- 
sions. Using the same Starburst99 models as above, we 
derive that the stellar contribution at 8/im ranges from 
1% to 30% with an average of 10% and the contribution 
at 24/im is negligible, less than 1% on average. 

The 8/im measurements compared to the short and 
long wavelength bands are summarized in Tabled where 
the ratio to the 4.5/im band is referred to as R\ and the 
ratio to the 24/im band as i?2- Increasing levels of PAH 
emission will result in decreasing values of R± and in- 
creasing values of i?2- We note that any correction for 
dust in the 3.6/im band will tend to increase the cal- 
culated 4.5/im dust emission and therefore increase the 
value of Ri for those galaxies without strong PAH emis- 
sion. 

The color measurements are presented graphically in 
Figure^ The figure shows that galaxies with high values 
of Ri tend to have low values of R2 , while those galaxies 
with a value of Rj_jn _the range of 0.01 to 0.09 (similar 
to our galaxy; cf. ILu 112004^ tend to have larger values of 
i?2, although with a large dispersion. We find that the 
separation of galaxies in this plot persists if we change 
our assumptions about how to compute the 4.5/im dust 
continuum, e.g., by chang ing the v a lue of a to 0.47 to 
match the value used by ILu et alJ l|2003fl or by using 
near-infrared data (such as the J band, which should 
not have a contribution from hot dust) to determine the 
stellar contribution to the 4.5/im band. This trend is 
consistent with the spectroscopic results summarized in 
the "PAH" column of Table ^ which demonstrate that 
all the galaxies with spectra free of PAH emission cluster 
in the upper-left portion of the diagram. 

To interpret thes e results, we used the SED models of 
iDale et aLl l)2001a|) . To compute photometry from the 
models (which do not contain a stellar emission com- 
ponent), we convolved them with the same spectral re- 
sponse curves discussed in § [21 The various flux density 
ratios were then computed as for the data. The emis- 
sion in the IRAC bands derived from these models is 
dominated by a fixed PAH spectral template, and as a 
result the models have a nearly constant value of R\, 
around 0.025. The models cover a range of R2 values 
that depends on the temperature of the thermally puls- 

6 http:/ /www. stsci.edu/science/starburst99/ 
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Fig. 1. — MIR colors of galaxies with a range of metallicitics. 
The 8/im band is compared to both shorter (4.5/xm, with stellar 
emission subtracted as described in the text) and longer (24/^m) 
wavelengths. Galaxies with known PAH features are indicated by 
filled markers, while galaxies known to lack PAH features are in- 
dicated by open markers. Error bars include a 10% relative cali- 
bration uncertainty between the IRAC and MIPS instruments and 
a 7% uncertainty on a, in addition to photometric uncertainties. 
The dashed line indicates our chosen separation between PAH and 
non-PAH galaxies in this color space. The range of points covered 
by the lDale et al.l 12001a!) SED models is shown as a solid line. 
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Fig. 2. — Galaxy metallicity as a function of the 8 to 24/^m color. 
Galaxies with colors or spectra that indicate that they have 8fim 
PAH features are displayed as filled markers, while galaxies which 
lack the 8fim PAH feature are shown as open markers. Squares 
and circles denote measurements with and without spectroscopic 
confirmation, respectively. The error bars on the 8/24/im ratio are 
the same as in Figure while the error bars on the met allicity are 
typically 0.05 or less (cf. Kobulnickv & Skillman 1996) and were 
all assigned an uncertainty of 0.05. Note that this figure does not 
require IRAC data (some of the 8(im measurements were made by 
MSX) and thus contains more points than Figure HI 



ing grains, along the locus of points outside the region to 
the upper left occupied by the PAH-free galaxies. Thus, 
the weakness of the 8/zm emission relative to the other 
bands appears to be due to the weakness of the PAH 
emission relative to the emission from small dust grains. 

Figure[5]shows the 8-to-24/im col or vs. metallicity (the 
solar metallicity on this scale is 8.7: IAllende Prieto et al.l 
1200 1|) . We label the galaxies according to the lack or 
presence of 8[im PAH emission using the color separation 
from Figure n They separate very cleanly according to 
metallicity, with a narrow transition region around 1/3 to 
1/5 solar metallicity. One additional galaxy, N GC 1569, 
has b een shown to have weak PAH emission l)Lu et al.l 
2003) and has an O/H-based metallicity 25% of solar, 
consistent with this overall behavior. 

The boundary between mctallicities with and without 
PAH emission is surprisingly sharp, and it seems likely 
that future studies will show more scatter. For example, 
there can be significant variations in metallicity within 
these systems (e.g., a factor of three is indicated for 
SBS 0335-052: lHouck et alJ2004bl) . Our characterization 
of this parameter with a single value is an oversimplifi- 
cation. In addition, "metallicity" includes enrichment 
of the ISM with a broad variety of elements, produced 
in different ways. Oxygen is produced predominantly in 
very massive stars; use of other elements produced in 
other ways as metallicity indicators may yield different 
behavior. Nonetheless, there is a robust trend toward 
low PAH emission in galaxies with relatively unpolluted 
ISMs. 

Why do the PAH features become weaker in low- 
metallicity galaxies? We reject extinction as a possible 
cause, both because the effect would have to (countcrin- 
tuitively) become stronger in galaxies with lower heavy- 
element abundance (and presumably less dust) and be- 
cause the direct measurement of extinction in the low- 
metallicity galaxy SBS 0335-052 shows it to be optically 



thin in the MIR ijHouck et alJl2004bH . It is also un- 
likely that starburst activity (i.e., the intensity of the 
star formation) is a major contributor to the effect: the 
correlation coefficient between Ri and the 24^m surface 
brightness (an indicator of the specific star formation 
rate, which we estimated from the integrated flux den- 
sities and the visible extent of the galaxies in the sam- 
ple) is much lower than the correlation coefficient be- 
tween i?2 and the metallicity (0.26 vs. 0.79). Also, while 
the well-known starburst galaxies in our sample (e.g., 
NGC 253, NGC 7714), do tend to have lower R 2 val- 
ues than other galaxies at similar metallicity, none of 
them have as small an R2 value as the lowest metallic- 
ity galaxies. One possible reason is the destruction of 
the PAH carrier in the harsh radiation fiel d of a low- 
metallicity galaxy l)Plante k, Sauvaeel 120 02). Another 
possibility is that these galaxies are truly young and 
the PAH molecule s have simply not had time to form 
l)Dale et al.ll200 1b). For example, it is well known that 
the ISM of a galaxy must be significantly polluted with 
oxygen (produced copiously in very massive stars) be- 
fore significant amounts of carbon are produced (by stars 
of a few M ) — this effect may delay PAH production 
(E. Dwek, private communication). The low- metallicity 
galaxies may simply lack the carbo n-rich AGB stars re- 
quired to form the PAH molecules l)Latterlll991|) . 

In a ddition, we find th at the PAH features included 
in the iDale et all l)2001aj) models fail to account for the 
high values of R\ prevalent at low metallicity (i.e., below 
8.0). The emission by dust in the models below 10/im is 
several orders of magnitude lower than the emission from 
the PAH template (which does not vary), so the simple 
remedy to the models of weakening the PAH spectrum 
does not allow them to match the low-metallicity data. 
The value of R\ derived from the model depends only 
on the shape of the PAH template spectrum and does 
not match the low-metallicity galaxies shown in Figure^ 



4 



(i.e., the model Ri value plotted in Figure ^ does not 
shift as the PAH features are weakened). Some other 
source of emission must be dominating i n these source s, 
such as hot dust in the 4.5/zm band (cf. ILu et al-l feOOS). 

4. CONCLUSIONS 

We have combined new MIR measurements from the 
Spitzer Space Telescope with results from the literature 
to show that the 8/tm to 24/im flux density ratio in star- 
forming galaxies depends strongly on metallicity. The 
8/tm to 24/xm color changes markedly between 1/3 to 1/5 
solar metallicity — the mean ratio below 1/3 solar metal- 
licity is 0.08 ± 0.04, while the ratio at higher metallicity 
is 0.70 ± 0.53. A MIR color-color diagram which com- 
pares the 8/im band to dust continuum measurements at 
4.5 and 24/mi shows that the change in the ratio is pre- 
dominantly due to a decrease in the 8/tm emission. We 
interpret this result as a weakening of the PAH features 
at low metallicity, an interpretation consistent with the 
spectroscopic evidence available in the literature. 

We show that existing empirical SED models cannot 
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account for the color change merely by weakening the 
PAH component of the models. In addition to the depen- 
dence of PAH strength on metallicity, SED models must 
account for another source of emission (e.g., hot dust) to 
explain the IRAC observations of low-metallicity galaxy 
colors. 

The shift in 8 to 24/im ratio occurs at a fairly high 
metallicity and will thus affect the spectra of many galax- 
ies. Future SED models to interpret Spitzer measure- 
ments at high redshift must take the dependence of 
metallicity on redshift into account. 
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TABLE 1 
Galaxy Colors and Metallicities 
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4.8e-2 


0.13 


0.05 


8.0 




1 

± 


"1 9 




UM 448 


l.lc-2 


5.6e-l 


0.03 


0.16 


8.0 


yes 




12 


9 


UM 462 


2.4e-3 


l.lc-1 


0.15 


0.06 


8.0 






2 




Mrk 930 


2.3e-3 


1.7c-l 


0.06 


0.06 


8.1 


no 


1 
± 


"1 9 


9 


II Zw 40 


l.le-2 


1.5e+0 


0.06 


0.07 


8.1 


no 


1 

1 


9 
Z 


13 


NGC 4670 


1.5e-2 


2.1c-l 


0.03 


0.36 


8.2 






14 




NGC 1156 




4.5e-l 




0.19 


8.2 




1 r i 

XD 


±u 




NGC 0598 




5.0c+l 




0.50 


8.5 




X t , -LO 






NGC 3077 


2.76e-l 


1.5e+0 


0.02 


0.48 


8.6 




1 


20 




NGC 300 


8.8e-l 


2.3e+0 


0.09 


0.70 


8.6 




21 


19 




NGC 2537 


4.1e-2 


2.4e-l 


0.02 


0.62 


8.7 




1 


22 




NGC 7714 




2.4e+0 




0.14 


8.7 


yes 


23 


20 


23 


NGC 2782 


3.8e-2 


9.6e-l 


0.02 


0.36 


8.8 


yes 


1 


14 


24 


NGC 4194 


6.5e-2 


3.1c+0 


0.03 


0.25 


8.8 


yes 


1 


20 


25 


NGC 5457 




l.Oc+l 




0.68 


8.8 




26 


19 




NGC 3031 


6.5e+0 


4.4e+0 


0.03 


2.02 


8.8 




27,28 


19 


29 


He 2-10 


6.0e-2 


4.9e+0 


0.02 


0.14 


8.9 


yes 


1 


11 


24 


NGC 253 




1.4e+2 




0.28 


8.9 


yes 


17,30 


19 


31 


NGC 5055 




6.1c+0 




0.77 


9.0 




17,30 


19 




NGC 7331 


9.5e-l 


3.6e+0 


0.01 


1.36 


9.0 


yes 


32 


19 


33 


NGC 224 




9.4e+l 




1.70 


9.0 


yes 


17,30 


19 


34 


NGC 5236 




4.2e+l 




0.50 


9.1 


yes 


1,30 


19 


29 


NGC 2903 


1.2c-l 


2.2c+0 


0.02 


0.50 


9.3 


yes 


1 


11 


9 
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